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We propose a device consisting in an antidot periodically driven in time by a magnetic field as 
a fractional quantum Hall counterpart of the celebrated mesoscopic capacitor-based single electron 
source. We fully characterize the setup as an ideal emitter of individual quasiparticles and electrons 
into fractional quantum Hall edge channels of the Laughlin sequence. Our treatment relies on a 
master equation approach and identifies the optimal regime of operation for both types of sources. 

The quasiparticle/quasihole emission regime involves in practice only two charge states of the an¬ 
tidot, allowing for an analytic treatment. We show the precise quantization of the emitted charge, 
we determine its optimal working regime, and we compute the phase noise/shot noise crossover as 
a function of the escape time from the emitter. The emission of electrons, which calls for a larger 
amplitude of the drive, requires a full numerical treatment of the master equations as more quasi¬ 
particle charge states are involved. Nevertheless, in this case the emission of one electron charge 
followed by one hole per period can also be achieved, and the overall shape of the noise spectrum 
is similar to that of the quasiparticle source, but the presence of additional quasiparticle processes 
enhances the noise amplitude. 

PACS numbers: 73.23.-b, 72.70.-|-m, 73.43.-f 


I. INTRODUCTION 

The on-demand single electron source (SES) based on 
a driven mesoscopic capacitor,^ has allowed to achieve 
interferometric experiments with individual electron and 
hole wave-packets propagating ballistically along integer 
quantum Hall (IQH) channels and opened the way to 
electron quantum optics.^ It relies on a quantum dot both 
tunnel-coupled to a quantum Hall edge channel and ca- 
pacitively coupled to a periodically modulated gate. In 
its optimal regime of operation, this periodic source emits 
exactly one electron in the first half-period and one hole 
in the second half-period.^^ 

This is achieved for an intermediate transparency of 
the point contact connecting the dot and the edge, when 
the gate is biased with a square voltage whose ampli¬ 
tude is equal to the dot level spacing. A complete 
characterization of this SES can be given in terms of 
a non-interacting picture in which the electron-electron 
Coulomb interaction on the dot can be effectively taken 
into account in terms of a renormalization of the level 
spacing. This allows to model the action of the meso¬ 
scopic CMacitor through the Floquet scattering matrix 
theory.™!] The SES allowed to perform Hanbury Brown 
and Twis^ (HBT) experiments with single electrons, as 
well as Hong-Ou-MandeP (HOM) collisions between two 
electrons propagating on opposite edge channels of the 
quantum Hall effect. 

While it is becoming clearer that electron-electron in¬ 
teractions in quantum Hall interferometric devices dra¬ 
matically affect the nature of the electronic excitations 
leading to their fractionalization,li2HIl] it jg natural to 
contemplate the fascinating possibility of operating such 
electronic interferometers in the fractional quantum Hall 



Figure 1. Artistic view of a (strongly asymmetric) antidot 
embedded into an Hall fluid and coupled with edge channels 
through tunneling amplitudes tn and tn respectively {\tL\ S> 
I till. Dark area represents the Hall fluid while the bright 
ones are in correspondence of the edges of the Hall bar (black 
arrows) and of the antidot (with black circles indicating the 
its energy levels). Blue dots indicated quasiparticles. The 
antidot is also pierced by a time dependent magnetic flux 


(FQH) regime,^ where interactions are at their strongest 
level. There, one would not be dealing with electrons, 
but rather with emergent Laughlin quasiparticle^!^ that 
carry fractional charge, and obey fractional statistics 
(anyons).l!^ The detection of fractionally charged quasi¬ 
particles in the FQH effect has been theoretically pre¬ 
dicted in calculations of the DC shot noise characteris¬ 
tics in a tunneling geometry between two counterprop- 
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agating edge states.^ Experiments gave a confirmation 
of these predictions.^ Alternatively, capacitive measure¬ 
ments using antidot geometries seem to also point toward 
the detection of fractional charges.^i^ 

Unfortunately there is so far no proposal for emitting 
such single quasiparticles along fractional edge channels: 
the simple picture of the source as a dot coupled via tun¬ 
neling to the edge cannot be easily generalized to the case 
of quasiparticle (QP) and quasihole (QH) emission in the 
FQH regime. Indeed, for dot transparencies in the opti¬ 
mal regime of emission, the Hall fluid is predicted to be 
so much depleted that only electrons can tunnel between 
the dot and the edge of the Hall bar.l^^ Conversely, at 
higher transparencies, the dot and the edge are strongly 
coupled so that the output signal reduces to the response 
of the edge to an applied voltage and is therefore not 
quantized in general. Granted, in this case, a properly 
designed Lorentzian voltage pulse could provide an effi¬ 
cient way to realize individual QP injection,l2SH2llat least 
for the simplest fractional states belonging to the Laugh- 
lin sequence.^ The main purpose of this work is to char¬ 
acterize the operation of a single Laughlin QP emitter 
which bears analogies with the driven mesoscopic capac¬ 
itor source, functioning with a driven antidot which is 
embedded in the fractional quantum Hall fluid. 

An antidot consists of a region of total depletion real¬ 
ized in the Hall fluid (see Fig. [^, and it can be driven 
either by modulating the Hall magnetic field or with a 
back gate. As mentioned above, this geometry has been 
considered as a tool to extract information about both 
the charge and the statistics of the QP excitations.^ In 
particular, earlier theoretical works investigated the adi¬ 
abatic pumping of the tunneling amplitudes connecting 
the antidot with the edges of the Hall bar as a way to em it 
a perfectly quantized fractional charge per cycle .1221^ In 
the stationary regime this peculiar geometry also offers 
the possibility to obtain a p ersist ent current induced 
by the piercing magnetic fieldMoreover, noise and 
higher order current cumulants present features able to 
disentangle the universal effects associated with the fill¬ 
ing factor and the non-universal ones related to the action 
of the external environment More recently, the same 

geometry has been discussed in the case of two dimen¬ 
sional topological insulators,!^ where the asymmetry of 
the antidot configuration turns out to be crucial in order 
to induce a spin polarized current across the sample.!^ 
On the experimental side, the periodicity of the conduc¬ 
tance peaks in this geometry has been measured as a 
function of magnetic field and back-gate voltage. This 
allows to extract charge and exchan ge statistics of the 
emergent excitations of the Hall fluid,“2Il5] even if these 
results have been longly debated.!^ 

In this paper, we show that a periodically driven an¬ 
tidot, can either behave as a single-QP source (SQS) or 
as a SES in the Laughlin sequence of the FQH regime. 
The setup is shown in Fig. where, in order to ensure 
the injection of fractionally charged QPs into the left 
(L) edge channel of the Hall bar only, a strong asymme¬ 


try is assumed. Here, for working purposes, the periodic 
drive is obtained by means of a time-dependent modu¬ 
lation of the magnetic field piercing the sample (equiv¬ 
alent modulation of the antidot levels can be achieved 
with a periodically biased voltage gate). In the present 
context, the oscillations imposed by the drive constitute 
only small perturbations with respect to the Hall quan¬ 
tizing magnetic field in order not to deviate from the 
fractional plateaus. The analysis of this system is car¬ 
ried out using the master equation approach,!^ however, 
in a different manner from what is usually discussed in 
the literature. Indeed, we need to face two new relevant 
problems: on the one hand the non-adiabatic time de¬ 
pendence of the drive has to be taken into account, and 
on the other hand we need to properly characterize the 
transient regime which is the physically relevant one for 
our purpose. 

By tuning the amplitude of the magnetic field oscilla¬ 
tions it is possible to induce the emission into the edge 
channel of a periodic train either of QPs and QHs, or 
of electrons and holes. The first case only involves two 
charge states of the antidot and can be solved analyti¬ 
cally. In particular, it is possible to compute the emitted 
charge per half-period as well as the associated current 
fluctuations, both these quantities showing remarkable 
resemblance with the i nteger q uantum Hall (IQH) meso¬ 
scopic capacitor setup.l^SEaiiU second case requires 
a full numerical treatment, due to the large number of 
charge states involved for the antidot. There, we observe 
a properly quantized electron charge per half-period, and 
a vanishing noise at zero frequency. 

However, the noise calculated at the drive frequency 
exceeds what would be expected for a SES, a conse¬ 
quence of spurious charge emissions randomly occurring 
during the half-period. All these tunneling processes do 
not affect the average charge or its fluctuations at zero 
frequency, but are inherent to the emission process of 
several quasiparticles in this system. Despite the great 
number of processes involved, in the regime of electron 
emission it is also possible to define an effective escape 
time for the excitations which can be experimentally ex¬ 
tracted from the measurement of the first harmonic of 
the current. 

The paper is divided as follows. In Sec. |H]we discuss 
the model of a finite length chiral Luttinger liquid de¬ 
scribing the antidot in the Laughlin regime of the FQH 
effect, coupled to the edges of the same fractional Hall 
bar. Sec. nn is devoted to the derivation of the tun¬ 
neling rates (for the master equation approach) as well 
as the relevant physical quantities, such as the occupa¬ 
tion number of the antidot, the current and noise along 
the edge, which are all essential in order to character¬ 
ize the performance of the device in both the SQS and 
the SES regimes. A discussion of the parameters for the 
optimal emission regime in the two cases is then car¬ 
ried out. In Sec. |IV| we consider the optimal regime of 
emission of QPs, where a simple analytical treatment is 
possible, characterizing the device as a perfect SQS. We 
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investigate the emitted charge, as well as the current fluc¬ 
tuations, showing that the considered setup behaves in 
exactly the same way as the SES realized in the integer 
regime, up to the renormalization of the charge of the 
carriers to account for their fractional nature, e* = ve. 
The optimal regime of emission of electrons and holes 
(SES) is investigated in Sec. using a full numerical 
approach. A perfect quantization of the average emitted 
charge 

is obtained also in this case, even though additional 
tunneling processes are required in order to properly de¬ 
scribe the system. These can be detected in the finite 
frequency noise and lead to additional complications in 
defining the escape time of the electron (hole) from the 
antidot. In Sec. |vg we perform estimates of the various 
physical parameters appearing in our SQS calculations 
in order to compare them with their counterpart in ac¬ 
tual experimental realizations in the IQH regime, and we 
justify the feasibility of our proposed setup. Sec. m is 
devoted to conclusions, while an Appendix discusses the 
connection between zero noise contribution at zero fre¬ 
quency and the absence of charge fluctuations during a 
period. 


II. MODEL 

The starting point of our discussion is the antidot ge¬ 
ometry proposed in Refs. [53 and [551 and schematically 
illustrated in Fig.[^ The FQH fluid is depleted by means 
of an electrostatic gate creating a circular “empty” region 
into the Hall bar. According to Wen’s hydrodynamical 
approactP^ for the description of FQH states belonging to 
the Laughlin sequenc^^^ with filling factor v = \/{2n-\-l) 
(n G N), the Hamiltonian associated with the boundaries 
(right edge i?, left edge L and antidot ad) is quadratic 
in terms of the edge-magnetoplasmon creation (aj and 
annihilation (ai^g) operators and can be written in the 
form (h = 1) 

+ O0 

= ( 1 ) 

S = 1 

where I = R,L, ad labels the various elements of the de¬ 
vice, e = 2TTvjLi is the energy associated with the plas- 
monic modes and E), = ttvv/Li is the charging energy 
associated with the zero modes, Li being the length of 
the Z-th edge. The propagation velocity v along the edges 
is assumed constant throughout all the setup. A) is the 
number of QPs enclosed by the edge I (with respect to 
a fixed background and iV$ = <i)/$o is the num¬ 
ber of elementary flux quanta d* of magnetic field pierc¬ 
ing the antidot section ($o = 27r/|e| the elementary flux 
quantum). Notice that the last term in Eq. Q is rem¬ 
iniscent of the minimal Aharonov-Bohm coupling j • A, 
with j the current density along the edge of the antidot 
and |A| = ^/L the vector potential felt by the anti¬ 
dot itself.^ In the following, we consider a finite length 


Lad = L only for the antidot with a consequent non-zero 
energy = Ac, while we assume the thermodynamic 
limit L[i,Ll —t - 1-00 for the R and L edges (thus yield¬ 
ing E^, E^ -G 0). Under these conditions the zero mode 
contribution plays a relevant role only for the dynamics 
of the antidot, while the other edges are only described 
in terms of their plasmonic modes. 

As long as the Hall fluid is present between the antidot 
and the boundaries of the Hall bar (see Fig. [^, the dom¬ 
inant tunneling process involves single-QPs with charge 
e* = zze (e the electron charge).!^ Therefore the local 
tunneling Hamiltonian connecting the antidot with the 
edges of the Hall bar is given by 


H 


V 




with j = L, R. (2) 


Here, the tunneling amplitudes tj are related to the over¬ 
lap between the Laughlin wave-function^^ on the differ¬ 
ent edges and decay exponentially with their distance. 
Tunneling processes occur at points Xj (see Fig.jl]), whose 
precise location is n ot rele vant as long as the tunneling is 
assumed to be local.li^*^ The vertex operator associated 
with the annihilation of one single-QP can be written, 
according to the standard bosonized description,!^ as 

= (3) 

V27ra 

The bosonic field ^pi{x) appearing in the exponent can be 
naturally decomposed into the sum of a plasmonic (yjf) 
and a zero mode (y)°) contribution given respectively by 

Stt 

=-j-vNix-xi (4) 

with a a finite length cut-off. The operator xi satisfies 


[xi,Lli>]=i6i^i> (5) 

and, once exponentiated, plays the role of a Klein fac¬ 
tor, which is essential to provide the correct exchange 
statis tical pr operties between excitations from different 
edges.1221281481 y^otice that the last phase factor in Eq. ^ 
has been introduced in order to satisfy the boundary con¬ 
ditions 


^iix + Li) = ^iix)E^^’'^^ (6) 

counting the n umber of fractional excitations enclosed by 
the edge.l2Sl2ll 

III. TRANSPORT PROPERTIES 
A. Tunneling rates 

To lowest order in the tunneling Hamiltonian the trans¬ 
port properties depend on the tunneling rates which, by 
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exploiting the periodicity associated with the antidot and 
assuming the plasmon mode s fully relaxed to thermal 
equilibrium, can be written 


+00 

p= — oo 

with e defined in Sec. El This combines the standard 
expression for the tunneling r ate a t finite temperature 
and infinite length of the edge^^^M) 


(dropping the subscript j in Tj because of the strong 
asymmetry), where the transition energies are given by 




{N - N^f 


{N' - N^f 


( 11 ) 


Since we want to tune the number of QPs on the antidot, 
we allow modifications of TV^, with respect to a reference 
value,!^ chosen to be a half-integer, so that 

Nfj) = Nq + - + S{t). (12) 


LUn 




(2^)2 
X B 


/3aJ< 

27r 


V Pi V .Pi' 
2"^*27r’2 *27r 


p/3C/2 


( 8 ) 


with the correction accounting for the finite length of the 
antidoiP^ 


T{u)p\ 


0(p-HO+) (9) 


evaluated analytically in the very low temperature limit 
pe:S> I {P being the inverse temperature) and for e/wc 
1. Here, we introduced the energy cutoff Wc = u/a, 
and used B[a,b] and r(a) as Euler’s Beta and Gamma 
functions respectively. These tunneling rates present a 
peaked structure near E/~ 0, with a maximum value 
decaying exponentially for negative energy (with a scale 
set by temperature), and as a power-law for positive one. 
Similar peaks also appear at EjE^ ~ {p G N°) as 

a consequence of plasmonic excitations, only with a less 
pronounced amplitude due to the damping factor Wp. 

When the tunneling rates are small compared to both 
the temperature and the charging energy (T^ < P~^ < 
Ec), it is possib le to restrict the analysis to the sequen¬ 
tial regimewhere only single-QP tunneling processes 
involving one excitation (incoming or outgoing with re¬ 
spect to the antidot) contribute to the dynamics of the 
system. In this approximation the relevant rates are the 
transition probabilities between the initial antidot state 
with N excitations at time t = 0 and the final one with 

± 1 excitations at time t. 


B. Master equation 


A simple and useful way to generalize this idea involves 
the master equation approach, which allows to character¬ 
ize the time e voluti on of the probability of occupation of 
the antidot.l^SHllIlZ] particular, the probability 7^7v(t) 
of having a fixed number N of QPs in a strongly asym¬ 
metric antidot (|ti| ^ |t/j|) at a given time t satisfies the 
first order differential equation 


dVN 

dt 




, ( 10 ) 


Setting n = N — No, p = N' — No and introducing Vn = 
VNljY^n+No^ the master equation can be conveniently 
rewritten under a matrix form as 


dVjt) 

dt 


= rvit), 


(13) 


where V{t) is a column vector whose elements are the 
occupation probabilities Vn(t), and P is a square matrix 
whose elements are given by 


fnp =r {-{n - p){n + p - l-26{t))Ec) 


{{n - q){n + q 

q& 


1 - 25{t))EP) 
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np' 


(14) 


In practice, while various possibilities are foreseeable, we 
will mostly focus on a square drive, defined over one drive 
period T as S{t) = S x Sgn — t). 


C. Occupation number, current and noise 

The occupation JV{t) of the antidot is readily obtained 
upon summing up the occupation probabilities solution 
of the master equation weighted by their corresponding 
number of QPs 


Af{t) = ^ nVn(t). (15) 

n^'L 


Of course, this occupation is defined with respect to the 
background reference set by No (cf. Eq. (Hi)- 

The total charge on the antidot is readily obtained 
from the occupation as Q{t) = e*N{t). It follows that 
the current flowing from the antidot to the edge is simply 
given by 


m = -^=-e*Y.Y.^t^pVp{t). (16) 

While the current provides crucial information on the 
operation of the antidot source, relevant information can 
also be obtained from its noise characteristics. Indeed, 
the finite frequency signal of the current current corre¬ 
lations allows a finer characterization of the operating 
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conditions of the QP or electron emitter, es pecial ly with 
regard to its dependence on the escape timeP^®^ 

When the escape time is much smaller than the period 
of the drive, the antidot emits QP or QH (alternatively, 
electrons and holes) in essentially a periodic manner and 
the main contribution to the noise is due to the uncer¬ 
tainty of the emission time within each period. This con¬ 
stitutes the regime of so-called “phase noise”, which is 
due to the random jitter of triggering of the drive and 
the actual emission time. On the opposite, when the es¬ 
cape time of QP (QH) (alternatively, electrons or holes) 
is much larger than the period of the drive, nothing guar¬ 
antees that the escape from the antidot really occurs, and 
these rare events give rise to a shot noise like contribu¬ 
tion. 

The master equation which is employed here allows to 
describe both regimes as well as the crossover between 
the two. However, prior noise experiments performed 
with the mesoscopic capacitor as the emitter tend (for 
convenience) to measure the noise at a frequency close to 
that of the drive, and chose instead to modify the escape 
time by tuning the transmission of the capacitor coupled 
to the edge, in order to explore the range of parameters. 

We define the current-current correlation as: 

Ciit,t') = {61{t)61it + t')), (17) 


with 


{Q{t)Q{t + t')) = e*'^ ^ nin2V{ni,t)'P{n2,t + t'\ni,t). 

ni ,712 

( 20 ) 

Numerically, this conditional probability is obtained by 
propagating the condition P{ni,t) = 1 through time, 
using the master equation in Eq. (13). 

Performing the same calculation for negative values of 
t', and accounting for the average charge, we have 


CQit,t') = {6Q{t)6Q{t + t')) 

= e*^ ^ nin 2 ^ 0 {t')V{ni,t)AV{n 2 ,t + t'\ni,t) 

ni ,712 

-I- 9(-t')V{n2,t + t')AV{nx,t\n 2 ,t + t')^ (21) 

with A.V{ni,ti\n 2 ,t 2 ) = (P(«i,ti|n 2 ,^ 2 ) - 'P{ni,ti). 
Computing the average over t and taking the Fourier 
transform, this ultimately leads to the following expres¬ 
sion for the frequency-dependent current noise 

i-oo pT 1, 

Si(uj) = 4e*^uj‘^ / dt'cos {ujt') / — 

Jo Jo ^ 

X nin2'P{ni,t)AV{n2,t + (22) 

Til ,712 


where SI{t) = J(t) — 

Because of the periodic drive, the current-current 
correlation depends on both t and t', and is 

T—periodic in time t. Since we are only interested in 
the behavior with respect to the time difference t', we 
consider a time-averaged quantity defined as 

5/(0 =2 / ^{6I{t)dI{t + t')) = 2{6I{t)6I{t + f)), 

Jo ^ 

(18) 

and the corresponding quantity in frequency space 

Si{oj) = J . (19) 

As it turns out, the simplest way to derive the cur¬ 
rent noise is to first access the charge noise as the two 
are trivially related: 5/(w) = a;^5(5(w).l^ Let us first 
focus on the charge correlation function {Q{t)Q{t -|- t')) 
and consider for simplicity that t' > 0. The charge cor¬ 
relation is only finite when the dot is occupied both at 
time t and at time t + t'. The actual value of the charge 
correlation is obtained by summing over all possible oc¬ 
cupations ni and 7 / 2 , multiplied by the joint probability 
of having ni QPs at time t and 7/2 QPs at time t + 

The latter is further written, using conditional probabil¬ 
ities, as the product of the probability V{ni,t) of having 
7/1 QPs occupying the antidot at time t, and the condi¬ 
tional probability V{n 2 , t + t'\ni,t) of having 7/2 QPs at 
time given that there were ni at time t. One is left 


D. Tuning the drive 

We now want to determine the optimal regime of op¬ 
eration of the source. First, we need to ensure that the 
antidot has sufficient time to emit/absorb QPs every half¬ 
period of the drive. To meet that goal, let us consider the 
situation of a rather large tunneling rate, say 70 T — 100 
where 



is the maximum tunneling rate between any two energy 
levels, and follow the evolution of the dot occupation over 
time for different values of the drive amplitude S. The 
results for the occupation ^(t) are presented in Fig. 
at filling factor 1 / = 1/3, for a given set of parameters 
satisfying the constraints of our model, namely PEf. ^ 1 
and EcfuJc ^ 1. These two plots clearly represent the 
two regimes of interest for the source. The upper panel 
shows variations of the antidot occupation between 0 and 

1, corresponding to the emission/absorption of a single 
QP. The lower panel shows variations between —1 and 

2, corresponding to the emission/absorption of an elec¬ 
tron charge (in the form of 3 QPs, consistently with the 
standard picture for the fractional states in the Laughlin 
sequence). Notice that the negative occupation is an arti¬ 
fact of the choice of the background reference, as the two 
extrema of the occupation are symmetric with respect to 

— Nq = 1/2. The optimal regime of operation cor¬ 
responds to a situation where the occupation is close to 
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an integer value at the end of every half-period, while 
behaving monotonously over each half-period. From the 
results of Fig. it seems the most appropriate choice for 
the drive amplitude is d = — 0.2 and 6 = — 1.2 for the 
SQS and the SES regime respectively. 

Strikingly, the evolution of the two types of occupation 
in the upper and lower panel (QP and electrons) bear 
strong similarities. When the drive amplitude is larger 
than the optimal drive in both cases, one sees that upon 
imposing the AC drive the occupation of the antidot over¬ 
shoots the zero or single QP or electron occupation: for 
QP (electrons), the trigger of the drive brings the antidot 
in a configuration which is slightly less than zero, while at 
the half period, the occupation is on average larger than 
one. This overshoot is slightly higher in the electron case 
than in the QP case. In the optimal case {6 = —0.2 and 
S = — 1.2 respectively), the response to the drive bears 
strong similarities with that of the mesoscopic capacitor. 




Figure 2. Occupation of the antidot over time as a function of 
the drive amplitude 5 for the SQS (top) and the SES (bottom) 
regime. Other parameters are: i/ = 1/3, Ec/lOc = 0.01, j3Ec = 
20 and yoT = 100. 


The remainder of this paper is devoted to the detailed 
study of these two specific cases. There, we will focus 
on the ’’stationary” state, still submitted to a periodic 
drive but independent of the initial conditions, namely 
the state of the system when the drive was turned on. 


IV. ANTIDOT AS A SQS 


We solve numerically the master equation in Eq. (10) 


over several drive periods, considering a vector column V 
of 31 elements (n = —15,..., 15) using a matrix version of 
the fourth-order Runge-Kutta scheme. As it turns out, 
most occupation probabilities are vanishingly small in 
this case, and one can focus on a very much reduced set 
of equations in order to properly describe the behavior 
of the system, allowing for an analytic treatment. 

In practice, one thus only needs to keep track of the 
probabilities for having a singly occupied or an empty 
antidot, namely Vi{t) and Voit) respectively. Focusing 
on a single half-period, corresponding e.g. to the emission 
of a QP, their dynamics is captured by the following set 
of equations 


dVijt) 

dt 

dVoit) 

dt 


= -I'Piit), 


(24) 

(25) 


where 7 = —fii = T{—2EcS). The other half-period, 
corresponding to the absorption of a QP, is described by 
similar equations, only exchanging Vq and Vi- Note that 
by definition, ^{6 = 0 ) = 70 and that 7 then rapidly 
decreases as one raises |(5|. 

This set of equations is simple enough to be propagated 
analytically. The occupation probability Vi{t) is thus 
given by 


P,{t) =Ai0)e 


- 7 * 


T 

for 0 < t < — 
- - 2 


A (i) = 1 + (Vi ( 0 ) - for ^ t < T. 

(26) 

while Voit) = 1 — Vi{t) at all times. The stationary state 
requires that A ( 0 ) = AA), so that Pi ( 0 ) = 2 coshlyr/i) ■ 
Numerical results for the evolution of these two proba¬ 
bilities over time are presented in Fig. [^for a given set 
of parameters at filling factor ly = 1/3. 


A. Dot occupation and current 


The occupation of the dot is obtained directly from the 
occupation probabilities and Eq. (151. Here this trivially 
reduces to the much simpler form 


j^it) = Ait), 

and similarly for the current 


lit) = -e 


dVijt) 

dt 


(27) 


(28) 


The evolution over time of these two quantities, evaluated 
numerically for different values of 70 T, is presented in 
Eig.g 
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Figure 3. Occupation probability of the two relevant levels of 
the antidot {Po and Vi) for the optimal regime of the SQS 
{5 = —0.2). Other parameters are: v = 1/3, EcjuJc = 0.01, 
j3Ec = 20 and 70 T = 100. The relevant scale entering the 
master equation is then 7 ~ 17.4512 T~^. 


Note that, because of the simple correspondence be¬ 
tween occupation of the antidot and probability of single 
occupation, the escape time r defined as the typical time 
associated with the emission/absorption process also cor¬ 
responds to the typical time-scale governing the evolution 
of the two relevant occupation probabilities. In other 
words, one trivially has r = 7 “^. 

The solution, Eq. (26), also allows to estimate the av¬ 


erage charge transferred to the edge during every half¬ 
period (in absolute value) 


Q = e* |A7(T/2)-A7(0)| 
= e* tanh ( • 


(29) 


confirming that the ideal operating regime is the one 
where 7 T 1 (Q —>• e*). 



Time t (in units of T) 


Figure 4. Occupation of the antidot (top) and current flowing 
between the source and the edge in units of ^ (bottom), in 
the optimal regime of SQS (5 = —0.2) for different values 
of 7 oT. Other parameters are: = 1/3, Ec/lOc = 0.01 and 
PE^ = 20 . 


which is very reminiscent of what was obtained for the 
SES in the integer quantum Hall regime 


B. Charge fluctuations 


The time-averaged charge fluctuations can readily be 
derived from the computation of the conditional occupa¬ 


tion probabilities following Eq. (21). 

The results are presented in Fig . | 5[ The y show that the 
time-averaged charge correlation CQ{t,t') vanishes expo¬ 
nentially with the same characteristic time scale as the 
dot occupation, i.e. 7 “^. The value taken at t' = 0 can 
be readily estimated from the expression for Vi{t) and 
reads 


CQ(t, 0) = e 


,2tanh(7r/4) 

7 T 


(30) 


One can verify that the charge fluctuations are thus given 
by the following exponentially decaying form 


7 T 


(31) 


C. Noise at the drive frequency 


Experimentally, the most accessible fluctuation-related 
quantity is the current noise probed at the frequency of 
the drive.l^HEI] \Ye computed this frequency-dependent 
current noise Si{uj) and evaluated it at the drive fre¬ 
quency Q = ^ for different values of JqT, or equivalently 
different values of the escape time r = 7 “^. The results 
are provided in Fig. 

Again, it makes sense here to compare the obtained be¬ 
havior to the one expected for the SES in the IQH regime, 
as the two are governed by a simila r set of equations in 
the present regime of operationl ^^ l ^^l Not surprisingly, the 
current noise obtained here reproduces exactly the result 
of the SES upon changing e into e*, namely 


5/(Q) 


~Y~ 


tanh 



1 + (Hr)" ■ 


(32) 


This excellent agreement means that there are no spuri- 



















































0.25 



Figure 5. Time-averaged charge correlation CQ{t,t') in units 
of e*^, for the optimal regime of the SQS (S = —0.2) and 
different values of joT. Other parameters are: ly = 1/3, 
Ec/uJc = 0.01 and f5Ec = 20. 



Figure 6. Frequency-dependent current noise Sj in units of 
* 2 

computed at the frequency of the external drive as a 
function of the escape time r from the antidot, for the optimal 
regime of the SQS {S = —0.2). The result is compared to 
the analytical expectation transposed from the known results 
of the SES in the IQH case, assuming a fractional charge 
e* = ve. Other parameters are: u = 1/3, EcjoJc = 0.01 and 
I3E^ = 20 . 


ous processes in the system, such as missed emissions 
compensated by double emissions, as well as QP/QH 
pairs. The source of the observed noise for r/T < 1 
is the incertitude in the time of emission of the excita¬ 
tions during the half-period, what is usually referred to 
as phase or jitter noise.l^ The considered emission pro¬ 
cess is therefore noiseless at long time-scales (see App. 

0 . 


V. THE ANTIDOT AS A SES 

We now turn to the case of emission/absorption of a 
single electron, by increasing the amplitude of the drive 
compared to the previous case of the SQS. Our results 


still rely on the numerical solution of the master equation 
over several drive periods, considering a vector column V 
of 31 elements (n = —15...15) and using a matrix version 
of the fourth-order Runge-Kutta scheme. 

Interestingly, this regime is much more complicated 
than the SQS. While one could have hoped to deal with 
only two occupation probabilities just as before (say V 2 
and V-i), it turns out not to be sufficient to properly 
describe the SES, which here cannot be reduced to a 
simple analytic treatment. 

In practice, even in the best case scenario (very large 
value of 7 oT) one needs to keep track of 8 different occu¬ 
pation probabilities in order to account for the behavior 
of the source over one drive period, and 6 if one focuses on 
only a half-period (the other 2 being recovered by sym¬ 
metry). Numerical results for the time evolution of these 
8 probabilities are presented in Fig. for 70T = 100. 
Focusing on a single half-period, corresponding e.g. to 
the emission of an electron, their dynamics is captured 
by the following set of equations 


—— = J^p with 


^ f _3__3 

0 0 

0 

0 

r-3,2 

f-2,-3 

r- 2,-2 0 

0 

r-2.1 

r- 2,2 

f-l.-3 

f_i._2 0 

f-1.0 

f- 1.1 

r- 1,2 

To.-3 

ho,-2 0 

fo,o 

fo,l 

ho ,2 

fl.-3 

0 0 

0 

fl,l 

ri.2 

V 0 

0 0 

0 

0 

^2,2 j 


(33) 


(34) 


where we omitted the time dependence for notational 
convenience. Here, the matrix elements T^p were defined 
in Eq. (14), and we kept only the non-vanishing con¬ 
tributions. The other half-period, corresponding to the 
emission of an electron, is readily obtained from the same 
set of equations upon exchanging P-3 o P3, P-2 ^ V4, 
P-i -O' P2 and Po -O' Pi (corresponding to the exchange 
between similar line styles in the top and bottom panels 
of Fig. [^. 

A few comments are in order from the observation of 
this reduced set of equations: 


• the apparent scattering that results in having the 
occupation of the antidot oscillating between 2 and 
— 1 does not involve the direct exchange of 3 quasi¬ 
particles, but rather processes involving between 1 
and 5 QPs, 

• during the emission process, ^2(1) follows an expo¬ 
nential decay, which introduces the characteristic 
time-scale T 2 = ^' 22 ^ 

• the above set of equations can be solved by succes¬ 

sively substituting the solution of the known prob¬ 
abilities into the next, according to the following 
order: ^2 — y 'P—n —^ 'Pi —^ 'P —2 —^ 'Po —^ 'P-i- 

This means in particular that all relevant occupa¬ 
tion probabilities can be viewed as a sum of expo- 
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Time t (in units of T) 




Time t (in xmits of T) 



Time t (in units of T) 


Figure 7. Occupation probability of the relevant levels for 
the optimal regime of the SES (<5 = —1.2), with 70 T = 100. 
Other parameters are: v = 1/3, EduJc = 0.01 and f5Ec = 20. 


nential terms of the form exp [note that all 

Tn^n are negative according to Eq. (14|]. 


This picture gets even more involved as one reduces 
the value of 70T. For 70T = 0(1), about 10 occupation 
probabilities are required in order to properly describe 
the electron source, signaling the importance of processes 
which involve transferring up to 9 QPs. 


A. Dot occupation and current 


The occupation of the dot is obtained directly from the 


occupation probabilities and Eq. (15), and similarly for 


the current. The evolution over time of these two quan¬ 
tities, evaluated for different values of 70T, is presented 
in Fig. 1^ Like the individual occupation probabilities, 
the antidot occupation can be written as a weighted sum 
of exponential terms of the form exp(r„_„f), as a direct 
consequence of Eq. (15). 


Note that, in contrast with the SQS, there is no proper 
way of unambiguously defining the escape time from the 
antidot. Following Eqs. (33)-(34), one might be tempted 


to introduce the characteristic time-scale T 2 associated 
with the exponential decay of V 2 , the probability for sin¬ 


Figure 8 . Occupation of the antidot (top) and current flowing 
between the source and the edge in units of ^ (bottom), in 
the optimal regime of SES (5 = —1.2) for different values 
of 7 oT. Other parameters are: = 1/3, Ec/lOc = 0.01 and 
/3Ec = 20. 


gle electron occupation of the antidot. However, these 
equations are not sufficient in the regime of low ■JqT, and 
V 2 no longer follows quite the same simple exponential 
decay. Similarly, one might want to define the escape 
time from the antidot Tad directly from the antidot occu¬ 
pation, e.g. from its initial rapid decay. However, since 
the occupation is a weighted sum of different exponen¬ 
tial terms, one can only extract estimates from the early 
rapid decay. 

A more convenient way is to follow the experimental 
procedure and extract the escape time from the first har¬ 
monics of the current. Indeed, one can write the first har¬ 
monics In = ^ Sq dtI as In = \In\e-’'^ and an ex¬ 
ponential decay of the current leads back to tanc/ = Hr. 
Transposing this to the present quantities, one can define 
the escape time 

T dtM{t) cos (nt) 

TO = - T —% -. (35) 

2’’’ Jp dt Af{t) sin (Of) 

The average charge Q transfered to the edge during 
every half-period was computed numerically from the ob¬ 
tained solution to the master equation and is presented 
in Fig. It qualitatively behaves as in the case of the 






























































10 




Figure 9. Average charge Q transferee! during one half-period, 
in units of e*, as a function of the escape time (corresponding 
to different values of Tdot, T 2 and rn) and for the optimal 
regime of the SES {5 = —1.2). Other parameters are: v = 
1/3, Ecjuc = 0.01 and PEc = 20. 


SQS, i.e. it shows a plateau at an integer value (here 3 
corresponding to the fact that one electron is equal to 
a bunch of 3 QPs in the Laughlin regime at v = 1/3) 
for the lowest escape times then decreases rapidly, con¬ 
firming that the ideal operating regime is the one where 
T <C T. Note however that it cannot be written as simply 
as before, and the form etanh(T/(4T)) (with any of the 
three definitions of r) though providing the good quali¬ 
tative behavior, does not fit exactly the numerical data. 


B. Charge fluctuations 


Figure 11. Frequency-dependent current noise Si in units 
*2 

of computed at the frequency fl of the external drive 
as a function of the escape time r from the antidot, for the 
optimal regime of the SES (5 = —1.2). Other parameters 
are: v = 1/3, EejoJe = 0.01 and PE^ = 20. The results are 
compared to the analytical expectation transposed from the 
SQS regime, Eq. (321, multiplied by 9 to account for the fact 
that the emitted charge now corresponds to one electron (3 
QPs). 


charge correlation CQ{t,t') vanishes rapidly with a char¬ 
acteristic time-scale which seems to be similar in value 
to Tdot) T2 and tq, without being quite exactly equal to 
any of those. This can be understood from Eq. (21), as 


CQ{t,t') appears as a sum of various exponential contri¬ 
butions, themselves obtained by combining two different 
terms of the form exp(r„^„t). 


The time-averaged charge fluctuations are derived 
from the computation of the conditional occupation prob¬ 
abilities. The results are presented in Fig. 

As already observed in the SQS case, the time-averaged 



Figure 10. Time-averaged charge correlation CQ{t, t') in units 
of e*^, for the optimal regime of the SES (<5 = —1.2) and 
different values of qoT. Other parameters are: v = 1/3, 
E^jiiJc = 0.01 and PE^ = 20. 


C. Noise at the drive frequency 


We consider now the frequency-dependent current 
noise Si{uj) associated with the SES and evaluate it at 
the drive frequency Q = ^ for different values of TqT, or 
equivalently different values of the escape time t (evalu¬ 
ated in three separate ways: ruot) T2 and tq). The results 
are provided in Fig. El 

It makes sense here to compare the obtained behavior 
to the one observed in the case of the SQS, correcting 
for the increased charge of the carriers. Since the noise 
is sensible to the charge squared, we plot for compari¬ 
son the result obtained in Eq. ( [3^ for QP injection in 
the previous section, multiplied by 9 (to account for the 
fact that the source now emits/absorbs 3 QPs per half¬ 
period). While the curves have qualitatively the same 
overall shape, the current noise of the SES remains much 
larger than the one for the SQS source, even when ac¬ 
counting for the increased transfered charge. This con¬ 
stitutes a clear indication of the importance of tunneling 
processes involving multiple QPs in this regime. 

Indeed, as argued when deriving the set of equations 
for the occupation probabilities, Eqs (33)-(341, a proper 


account of the behavior of the SES requires to consider 



































11 


tunneling processes involving up to 9 QPs whereas we 
only emit an average of 3 QPs or less per half-period. 
Although these processes are responsible for the signif¬ 
icant increase of the current noise at finite frequency, 
thus contributing to the phase (finite frequency) noise, 
they do not affect the quantization of the emitted charge 
per half-period (see Appendix]^ so that the numerically 
evaluated frequency-dependent noise still vanishes at zero 
frequency. 


VI. ESTIMATING THE PHYSICAL 
PARAMETERS OF THE SQS 

For our simulations, we have chosen specific values for 
the various parameters which satisfy the constraints of 
our model, but may otherwise look random. Here we 
show through a simple estimate of the scales at play in 
an actual realization of our setup that these quantities 
are actually compatible with what is typically observed 
in electron quan tum op tics experiments carried out in the 
integer Hall caseIn what follows, we reintroduce the 
proper factors of h and fcs when needed. 


A. Energy scales 

Focusing on the state at filling factor i/ = 1/3, assum- 
ing a small antidot of circumference L « 1 /rrrP 2 and a 
propagation velocity along the edge v « 10 ^ m/s, one has 
the charging energy 


EjkB « 0.8K (36) 

(expressed in units of temperature) which is smaller, but 
of the same order of magnitude as the energy gap ob- 
ser ved in the SES in the IQH regime (A « 1.4 to 4.2 
K). bFF° l The plasmonic energy is larger 

e/ks ~ 4.8K, (37) 

but plays only a marginal role in the dynamics of the 
device. 

The temperature of the system, as assumed in the sim¬ 
ulations, is 


0 = 0.05EjkB « 40mK. (38) 

It is comparable to the one obtain ed in experiments in¬ 
volving the antidot geometrjEUlll and not out of reach 
from nowadays experimental techniques which allow to 
carry out el ectro n quantum optics measurements (0 « 60 
to 100 mK). 12112] 

The optimal operation regime of the SQS discussed in 
the text is reached for a drive amplitude S = — 0 . 2 , which 
corresponds to a magnetic field fluctuation of 

AB = 0.2 « lO-^T. (39) 


This represents a small and experimentally achievable 
fluctuation with resp ect to the magnetic field needed to 
realize the FQH state^^MU (typically of the order of lOT, 
depending on the two dimensional electron density of the 
considered sample). 

The high-frequency cutoff uJc = v/a represents the 
highest energy scale in the model. We assume that the 
short length cutoff a is typically set by the magnetic 
length of the system, so that a « 10 nm. This in turn 
leads to Wc ~ 10 THz, or equivalently 

hojJkB « 76K (40) 

from which one readily sees that Ec/oJc — 0.01 as used 
in our simulations. 


B. Time scales 


In order to operate the antidot in the optimal regime 
of emission, the typical time-scales of the device need to 
satisfy the relation 

T 

To <-f~^ < - (41) 

where tq = L/v is the time required to make a loop 
around the antidot. The first inequality ensures the va¬ 
lidity of the continuous limit assumed in the master equa¬ 
tion approach, as tq corresponds to the discretization 
time associated with the semiclassical processes.]^ The 
second one guarantees that the emission of a QP (QH) is 
achieved during one half-period. According to the previ¬ 
ous estimation one has 

L , , 

To = — « lOps. (42) 

V 

Assuming a driving frequency / = 500 MHz, close to 
the GHz regime investigated in Refs. [T]and|38l the drive 
period is given by 

T=j^2ns (43) 

so that the half-period is of the order of the nanosecond. 

This enforces the emission time from the SQS to vary in 
the range 10 to lOOOps, compatible with what is currently 
observed in similar experiments (y”^ ~ 60 to 900ps i n the 
optimal regime for the SES in the IQH case) Note 

that any value of the emission time 7 “^ in this range also 
satisfies the condition 


H'j < kBO < Ef. (44) 


necessary to _ Rlfill the sequential tunneling 

approximation.1221551 

The maximum value of the emitted current is directly 
related to the emission time, as one has 


, exp( 7 T/ 4 ) 

"^2cosh(7r/4) 


(45) 
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which corresponds to a value of the current between 
lOpA and InA, a reasonable range in comparison with 
what is experimentally measured for the SES in the in¬ 
teger regim^ and for the continuous current measure¬ 
ment in the weak backscattering regime at fractional fill¬ 
ing factor!^ 

Finally, using our energy-scale estimates, we can prop¬ 
agate these bounds on the emission time and obtain 

10"'^ < |tLp < 10■^ (46) 

for the range of transmission amplitude between the an¬ 
tidot and the nearby edge channel. 


VII. CONCLUSIONS 

In this paper we have analyzed a strongly asymmetric 
antidot geometry realized through depletion of the Hall 
fluid, and periodically driven in time by a modulated 
magnetic flux. Through a master equation approach, we 
discussed the possibility to use this kind of setup both 
as a SQS and a SES. In the first case, only two charge 
states of the antidot are involved and the dynamics of 
the system allows for a tractable analytic treatment. It 
shows the quantization of a fractional charge e* = ve 
emitted during each half-period and fluctuations analo¬ 
gous to the ones observed for the electron emission in the 
integer quantum Hall case. The electron emission regime 
proves more complicated as various charge states of the 
antidot are involved, requiring a full numerical treatment. 

Here, we observe that we can reach a regime where the 
emitted charge over a half-period is precisely e (with the 
precise opposite charge emitted over the second half of 
the period). However, the noise measured at the drive 
frequency is strongly enhanced with respect to what is 
observed in the integer regime due to the random emis¬ 
sions of additional excitations having zero mean charge 
(quasiparticle-quasihole pairs, etc.) which provide an es¬ 
sential contribution to the dynamics of the system at 
finite frequency. In spite of the presence of these ad¬ 
ditional tunneling events, it is possible to extract infor¬ 
mation about the escape time of the electron by looking 
at the first current harmonic as is usually carried out in 
experiments. 

It is worth underlying that in previous works, 
discussing similar geometries realized with Luttinger 
liquid^^ (typically in CNli^or two dimensional topolog¬ 
ical insulator^^ the role of the magnetic flux is played 
by an external gate voltage. The possibility to exploit a 
fluctuation of the same gate voltage used to realize the 
antidot also to induce the QP injection could be fruitful 
from the experimental point of view, an external elec¬ 
trostatic gate being easier to tune than a magnetic field. 
However, in this case, modifications of the antidot geom¬ 
etry are expected and their effects on the functionality 
of the presented device have to be carefully taken into 
account. 


The present study could prove quite relevant for the 
implementation of both theoretical and experimental in- 
vestirations of interferometric phenomena in HBT^^ and 
HOIVP setups in the context of electron quantum optics. 
Indeed, the present SQS or SES source could be embed¬ 
ded in a quantum Hall bar which is divided in two by 
a central quantum point contact, where either the HBT 
partitioning of the source could be analyzed in order to 
quantify the production of spurious excitations, or alter¬ 
natively with two antidot sources where two QP colli¬ 
sions could be achieved in order to probe the overlap of 
QP wave packets. In the near future, we are determined 
to model such experiments assuming that either ideal 
QP wavepackets or ideal electron wavepackets have been 
deposited on the fractional Hall edge,l^ but a proper de¬ 
scription of the source, such as presented in the context 
of this work, will ultimately be necessary to bring the 
description sufficiently close to experimental reality. 

Furthermore, the present study was achieved assuming 
the weak coupling hypothesis where tunneling rates are 
smaller than the electronic temperature. The numbers 
we provided in Sec. |VI A| seem to point out that this 
hypothesis is justified in the context of present experi¬ 
mental working conditions (fcs© ~ 50mK). However, 
upon either increasing the tunneling rates (by bringing 
the antidot closer to the edge, which would increase the 
QP tunneling amplitude) or by working at much lower 
temperatures, a coherent description of tunneling will be 
eventually required. This constitutes a truly challenging 
task, which may have to rely on Keldysh non equilibrium 
Green’s function formalism in order to describe the time 
dependence of the current, while taking into account the 
(finite) QP occupation of the antidot in its evolution. 
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Appendix A: Charge fluctuations during a period 

This Appendix discusses the relation between the cur¬ 
rent noise at zero frequency and the fluctuations of the 
charge emitted by a SQS or a SES during one period. 
Let us start by considering the general definition of the 




13 


noise 


SM = {I{t)I{t'))-{I{t)){Iit')) (Al) 


where / is the current flowing along the considered chan¬ 
nel and the averages are taken with respect to an arbi¬ 
trary initial state. 

The operator associated with the charge emitted dur¬ 
ing one period T is given bjisa 


Q{T) = 


r+T/2 


l-T/2 


drlij). 


(A2) 


the first term in Eq. (A31 becomes 
r+Tl2 


l-T/2 


dtdt'S(t, t') = 


r+T/2 

l-T/2 

r+T/2 

l-T/2 


dt 


dt 


+ 00 


' —OO 
/• + 00 


drS (^t + 


+00 

dT 

m— — oo 


T - 



1 

T 


X 



-iuirgim), ) 

27r ^ ’ 


= S^m=o\uj = Q), (A6) 


where in the second line we introduced the Fourier trans¬ 
form associated with the variable r and the series rela¬ 
tive to t, while in the last line the only remaining term 
is given by the m = 0 harmonic of the noise evaluated at 
zero frequency. 

Summarizing the above results, one finds 


Together with Eq. (Al), the above definition leads to 


(Q(T)Q(T)) = 



dtdt'{I{t)I{t')) 

dtdt'S{t,t') + {Q{T)){Q{T)). 

(A3) 


(Q(T)Q(r)) = (a; = 0) + {Q{T)){Q(T)) . (A7) 

Concerning the cases considered in this paper, analytical 
and numerical evidences suggest that 

^(m=0)(^ ^ 0) « 0 (A8) 

and consequently that 

{Q{T)Q{T))-{Q{T)){Q{T))^Q. (A9) 


Taking into account the periodicity of the noise signal 
and using the standard change of variable^ 


t — t+ - (A4) 

= (A5) 


This is a clear signature of the absence of charge fluc¬ 
tuations during one period for both the SQS and the 
SES realized in the fractional regime, in full analogy with 
what is observed for the SES in the integer quantum Hall 
case.l^ However, charge fluctuations due to the uncer¬ 
tainty in the moment of QP emission from the source 
still remain. This so-called phase or jitter noise provides 
an intrinsic finite frequency contributiorP which can be 
detected by measuring the noise at the frequency of the 
drive. 
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